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Abstract 

It is shown that the azimuthal dependence of the distribution of hadrons in 
a quark jet is a probe of the transverse spin of the quark initiating the jet. 
This results in a new spin-dependent fragmentation function that acts at the 
twist-2 level. One example of a process where it contributes is semi-inclusive 
deeply inelastic scattering with a transversely polarized nucleon target but with 
an unpolarized electron beam. This process is treated in detail. Another process is 
the cross section for singly polarized hadron-hadron scattering when two high p± 
hadrons are measured in the final state and are close to back-to-back in azimuth. 
The new fragmentation function is sensitive to the coupling of the fragmentation 
process to (spontaneous) chiral symmetry breaking. 

1 Introduction 

An important challenge to QCD theorists is to devise methods of measuring the polarization 
state of a parton coming out of a hard scattering. In the case of longitudinal polarization 
for the parton, Nachtmann [1] showed how a certain three-particle correlation within a jet 
could be used. He suggested several processes where it could be measured, in particular 
deep inelastic neutrino scattering. Later, Dalitz, Goldstein and Marshall [2] and Einhorn [3] 
showed how to probe the helicity of a heavy quark. 

Recently Efremov, Mankiewicz and Tornqvist [4] rediscovered the Nachtmann idea, which 
they called the 'handedness' of a jet, and they showed how to measure it in e + e~ annihilation. 
They also considered the possibility of probing the transverse polarization of quarks. This 
idea was independently discovered in [5]. 

In this paper, I will present another technique sensitive to the transverse polarization of a 
quark. One measures the single particle distribution as a function of transverse momentum, 
in a process such as semi-inclusive deeply inelastic scattering, where the jet axis can be 
precisely defined. The initial-state hadron is transversely polarized, and the hard scattering 
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provides a spin transfer to the final state. There should be a spin asymmetry in the azimuthal 
distribution of the outgoing hadron about the jet axis. 

One could argue that the large number of hadrons produced in a jet would completely 
dilute this asymmetry. However, it is known that for flavor quantum numbers, the leading 
particle in a jet is correlated in flavor with the quark initiating the jet [6]. Moreover, much 
of the hadron production results from the emission of gluons at small z; this preserves the 
helicity state of the quark, and it is the quark that is responsible for the leading hadrons. 

Measurements like this probe three quantities, none of which have yet been measured: 
the transverse spin of a quark in a transversely polarized hadron, the transverse spin depen- 
dence of hard scattering, and the spin-dependence of fragmentation. The difference between 
distributions for the transverse spin and for the helicity of a quark in a hadron is a measure 
of relativistic effects in the hadron's wave function [7]. Transverse spin dependence of the 
fragmentation involves a helicity flip for the quark, and is therefore sensitive not only to the 
spin dependence of fragmentation in general but to the presence of (non-perturbative) chiral 
symmetry breaking. 

For the single-particle measurement, one must have a process where there is a definite 
axis or plane with respect to which transverse momentum of an outgoing particle can be 
defined. Two obvious examples are: 

• Semi-inclusive deeply inelastic scattering, e + p 1 ' — > e + H + X , when the kinematic 
variables xbj and z are not small, and the transverse momentum of the observed hadron 
(relative to the plane of the initial proton and the two electrons) is moderate. This is 
the process that I will treat in detail in this paper. 

• p+p 1 * — > H1 + H2 + X, when the two observed final-state hadrons have large pt relative 
to the axis of the initial-state particles and are nearly back-to-back in azimuth about 
that axis. 

Here, p 1 ' denotes a transversely polarized hadron. Undoubtedly, other similar observables 
can be easily devised. It is also standard to define a jet axis by calorimetric measurements of 
the jet for example. In that case one can measure the azimuthal distribution of the leading 
pions about the axis. This probably suffers from lower precision in the definition of the axis. 

In the proton-proton case, if one makes measurements at sufficiently large center-of- 
mass rapidity, then one can have the collision of a small x gluon and a moderately large 
x valence quark. This allows one to take advantage of the large number of small x gluons 
and the expected large polarization of the valence quarks. This would give both a large spin 
asymmetry and a large cross section. 

Both the single-particle measurement suggested in this paper and the two-particle mea- 
surement suggested in [4,5] are twist-2: that is, the spin asymmetries are not suppressed by 
a power of the large momentum scale Q in the hard scattering. This is important because 
of the paucity of useful twist-2 asymmetries involving transverse spin. The problem is that 
transverse spin asymmetries involve off-diagonal elements in the helicity density matrices of 
the quarks, and that the hard scattering preserves quark helicity. (It is quarks, not gluons 
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that are relevant here [8].) Thus to get a twist-2 asymmetry, the non-perturbative part of 
the process (consisting of the parton distribution and fragmentation functions) must contain 
an even number of factors that flip quark helicity. 

The obvious measurements that are sensitive to transverse spin at the twist-2 level involve 
collisions of two polarized hadrons, and to get a large asymmetry one wants to be in the 
valence region for the quarks entering the hard scattering. The two standard processes 
are jet production and the Drell-Yan process. But jet production has a numerically small 
asymmetry [9] in the hard scattering. The Drell-Yan process, with a large asymmetry in the 
hard scattering [10], needs a polarized antiproton beam with good luminosity to make it a 
valence process with a reasonable event rate. 

In contrast, the spin asymmetries discussed in this paper, and those in [4,5] have their 
second measured spin in the final state instead of the initial state. It is interesting that these 
fragmentation measurements are easier for transverse polarization than for longitudinal. A 
transverse polarization can be probed by measuring two suitable momentum vectors (e.g., 
two particles in a jet, or one particle and a jet axis). But a longitudinal polarization requires 
three vectors, to define a handedness — see [1,4], and also the work of Carlitz and Willey [11]. 

The new spin dependent fragmentation effect that is the subject of this paper gives a 
leading-twist contribution only to processes that are sensitive to intrinsic transverse momen- 
tum. But it may also have consequences for other processes. For example, consider the 
inclusive cross section for production of a high pj_ hadron in hadron-hadron collisions. When 
only one of the initial hadrons is transversely polarized, QCD predicts that the resulting 
spin asymmetry is of twist-3. There is a notorious contrast between this prediction and the 
large asymmetry actually measured [12]. The experiments are admittedly at rather modest 
energies (fixed target with beams of up to 200 GeV). Methods for treating twist-3 effects are 
currently begin developed [7,13,14]. 

In that process, there is effectively an average over intrinsic transverse momentum, and 
in particular over its azimuthal angle about the jet axis. Therefore the spin-dependence of 
the fragmentation that is treated in the present paper gets canceled, at the leading twist 
level. But the cross section is a steeply falling function of the hadron's p±, so that the 
azimuthal average need not work very well. The result could be a numerically large twist-3 
contribution in the fragmentation. The authors of [14] have not yet treated twist 3 effects 
in fragmentation. 

The layout of this paper is as follows. I define the fragmentation and distribution func- 
tions in Sec. 2. Then in Sec. 3 I show how these get used in deeply inelastic scattering, and 
in Sec. 4 I explain the polarization sensitive measurements that can be made; that section 
is perhaps the one of most use to experimentalists. In Sec. 5, I make a few remarks about 
hadron-hadron scattering. Then I exhibit some simple model calculations in Sec. 6 to show 
that the fragmentation asymmetry I propose does not violate any fundamental principles. I 
summarize my conclusions in Sec. 7. 
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2 Spin-Dependence of Parton Distribution and Fragmentation Functions 

In this section, I will first review the formal definitions of the parton distribution and frag- 
mentation functions when there is a measured transverse momentum. Then I will show how 
to extend the definitions to treat nontrivial polarization. These quantities will get used in 
factorization formulae for the cross section, as explained in later sections. 

In the usual factorization theorems [15,16], one works with parton densities integrated 
over transverse momentum. But when one has a cross section with a measured small trans- 
verse momentum variable, one must use the unintegrated distributions. In QCD, there are 
some interesting effects associated with Sudakov form factors, that make the resulting fac- 
torization theorems quite nontrivial [17]. The Sudakov effects are spin-independent, and 
we will not bother making them explicit here, since our purpose is to examine the novel 
effects associated with polarization. However when the energy of the experiment increases, 
the Sudakov effects will dilute our asymmetries by smearing out the transverse momentum 
distributions. 

We will denote the unpolarized distribution and fragmentation functions by fi/n and Dh/% 
respectively, when transverse momentum is integrated over. To denote the corresponding 
quantities with unintegrated transverse momentum, we will use the same symbols, but with 
a hat over them: fi/ H and Dr/i- 

2.1 Parton Distribution Functions 

We define parton distribution functions by formulae motivated by light-front quantization. 
These quantities are precisely those that occur in the factorization theorems [18,19,20]. 

It follows from the parity and time-reversal invariance of QCD that the number density 
of quarks is independent of the spin state of the initial hadron, so that we have 

fa/A(x, \k ± \) = J d y~^y± e - ixp +y- +ik± . y± ^ (Q; y - y±) 1_ {p y (l) 

We have ignored here the subtleties needed to make this a gauge invariant definition: an 
appropriate path ordered exponential of the gluon field is needed [18]. The coordinate 
frame in which this definition is applied is one in which the hadron \p) has zero transverse 
momentum: p± = 0. 

Sivers [21] suggested that the k± distribution of the quark could have an azimuthal 
asymmetry when the initial hadron has transverse polarization. However, such an asymmetry 
is prohibited because QCD is time-reversal invariant. This is shown in the appendix. 

As explained in [10,22], we must consider the quark (or gluon) a to be equipped with a 
helicity density matrix. Since QCD is invariant under parity and time reversal, the density 
matrix for a quark differs from unity only if the initial hadron A is itself polarized. Then 
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the transverse spin asymmetry of a quark is defined by: 

s^A Ta/A f a/A (x, k±) 



where Sj_ is the transverse part of the spin vector of the initial hadron, normalized so that its 
maximum size is unity: < 1. In eq. (2), I have used the notation of [22], where definitions 
are given for the case that the transverse momentum is integrated over; the definitions given 
for that case given by Jaffe and Ji [7] differ only in notation. I have defined to mean the 
ratio of quark polarization to hadron polarization; it is a kind of asymmetry or spin transfer 
function, and in general will depend on x and k±. Then (with a subscript T) means the 
parton distribution weighted by the transverse spin asymmetry. 

Similar definitions can be given for the distribution of gluons. It would be interesting to 
work out the details. For the deep-inelastic process treated in this paper, we will see that 
we will not need the definitions for gluons. 

One can also write helicity asymmetries. But we will not need them in this paper, because 
we will work with fragmentation observables that are not sensitive to quark helicity. 

2.2 Fragmentation Functions 

Fragmentation functions with transverse momentum are defined in a similar fashion to 
the parton distribution functions. An important difference is that the observed final-state 
hadron, assumed to be a pion, has no polarization, so that a spin-transfer function analogous 
to eq. (2) does not exist. On the other hand we can no longer apply the generalization of 
the time-reversal argument that prevented an azimuthal asymmetry in the number density, 
as we will explain. 

The unpolarized fragmentation function, to find a hadron H in the decay products of a 
quark of flavor a is [18] 

D H/a (z,k L )^Y,J ^^e^"*^ (3) 

The factor of 1/12 is the product of a factor 1/2 which occurs with all these definitions for 
fermions and a factor 1/6 for an average over the spin and colors states of the quark. The 
sum over X is over all final states containing the chosen hadron. 

This definition is easily generalized to give the transverse spin dependence AD of the 
distribution of hadrons in a polarized quark of transverse spin Sj_: 

AD H/a (z, k±, s±) = 

(4) 

5 



slfTa/A^X, k±) 

dy-d 2 y ± _ 



I 
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^y-^-y- (p| ^ a (o, y - y± ) Va(0) \p), 



This is permitted to depend on the spin through the following scalar quantity: 

e KXfj,u s ±kj_Pjjn u , (5) 

where n M = 5t is the vector that is used to define the light-front momentum fraction z: 
z = pjj/k + = pn ■ n/k ■ n, with /c M being the momentum of the quark. (Note that the above 
definitions imply that k~ is integrated over, while k± is the transverse momentum of the 
quark relative to the hadron H .) 

There is some disagreement in the literature on the normalization of the fragmentation 
functions. The definitions given above are the most convenient for theoretical analyses since 
the bilocal vertex for the quark depends only on momenta defined at that vertex. However, 
a probability interpretation is most naturally given in terms of the transverse momentum 
of the hadron relative to the quark, which is p± = —zk±. Then the probability density in 
z and pj_ is z~ l D(z, —zpj_). As for the density integrated over transverse momentum, i.e., 
D(z) = / d 2 k±D(z, k±), the probability density is zD(z). (This implies that the momentum 
sum rule is J2h Jo z2 Dn/a( z ) dz — 1.) 

Now time reversal invariance prohibits an asymmetry proportional to (5), unless there 
are nontrivial phases from final state interactions; but such phases can surely exist in a 
strong interaction fragmentation. Another way of saying this comes from observing that a 
combination of time reversal and parity invariance transforms out-states to in-states of the 
same momentum content. (Compare the arguments given in the appendix for the distribution 
functions.) If the in-states were the same as the out-states, then we could show that the 
asymmetry proportional to (5) is zero. For strongly nonperturbative interactions, such as 
are involved in hadronization, this is surely not true. Similar arguments have been applied 
to polarized pion-nucleon scattering [23]. 



3 Semi-Inclusive Deeply Inelastic Scattering 

3.1 Kinematics 

Semi-inclusive deeply inelastic scattering is the process e + A^e' + B + X where A is some 
arbitrary initial hadronic state of momentum p A) and B is an observed hadron in the final 
state of momentum p^. The momentum transfer from the lepton system is q^, as usual, and 
we define the masses of hadrons A and B to be Ma and M B . 

In the absence of polarization, there are four Lorentz invariant variables for the hadronic 
system, viz., q 2 , Pa ■ q, Pb ■ q, and Pa - Pb- A useful set of variables was defined by Meng, 
Olness and Soper [24]. The two variables not involving the final-state hadron B are, as usual, 

and 

Q = \f Z Q 2 - (7) 
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Fig. 1. Kinematics of semi-inclusive deeply inelastic scattering. 
Then as in [25] one defines 

z = P -^. (8) 
Pa ■ q 

Finally, Meng et al. define as the transverse momentum of g M , in a frame where yf A and 
Pb have zero transverse momentum. Up to corrections that are irrelevant in the Bjorken 
limit, this means that 

<£ = <f + aup£--p£, (9) 

where x A = ~Pb • q/pa-Pb = x B - s (1 - \q\\/Q 2 )- 

The deep inelastic region is where Q is made large, with x-q] and z held fixed and not 
close to their endpoints and 1. We will always assume in this paper that the scattering 
is taken to lowest order in QED, with a single photon being exchanged between the lepton 
and the hadronic system, Fig. 1. 

The reason for defining the variables xbj, z and q± is that they have a simple interpretation 
in the parton model. There, it is assumed that the dominant contributions to the cross section 
have the form of Fig. 2. The virtual photon interacts in Born approximation with a single 
quark, which is close to its mass shell and which has low transverse momentum on the scale 
Q. Then when hadron B is part of the 'current quark jet' produced in the hard scattering, 
z has the interpretation of the fraction of the jet's momentum that is carried by the hadron. 
As usual, XBj has the interpretation of the fraction of the momentum of the incoming hadron 
A that is carried by the parton that enters the hard scattering. 

To treat intrinsic transverse momentum for the initial state and for the fragmentation, 
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Fig. 2. Parton model for semi-inclusive deeply inelastic scattering. 
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we need a suitable frame in which to define them. First we define what we will call the 
'parton model jet axis': 

^EE^ + XBj^. (10) 

This would be the jet momentum if there were no intrinsic transverse momentum. Even 
in the presence of transverse momentum, the definition (10) gives a convenient axis with 
respect to which to define the transverse momentum of the outgoing hadron B. 

As stated above, q± is the transverse momentum of the virtual photon in a frame in which 
the two hadrons A and B have zero transverse momentum. This is perhaps the simplest 
frame in which to derive the factorization theorem at low transverse momentum. But when 
we consider the azimuthal distribution of the outgoing hadron B it will be more intuitive to 
consider the transverse momentum of B relative to the parton model jet axis. This we can 
define by 

Pbt=Pb~-Pj- (11) 

Note that to avoid some possible confusion we use the subscript t here rather than j_. This 
is to emphasize that pg T is a four vector rather than the transverse two-component vector 
defined in a particular frame: the frame in which we define the fragmentation function will 
have zero transverse momentum for B. 

Up to corrections of relative order q"j_/Q 2 , the transverse momentum is —zq^_. 

3.2 Factorization without intrinsic transverse momentum 
In QCD, there is a standard factorization theorem: 

E ' EB dn>d*p B = ^S di Sj fa/A ^ ] (Hard scatterin g) D B/b(0, (12) 

which is correct up to power-law corrections. When the hard scattering factor is replaced by 
its lowest order approximation, we regain the parton model. In normalization the hard scat- 
tering function is the invariant cross section for lepton-parton scattering with the outgoing 
parton momentum set to Pb/C- 

This factorization theorem represents the asymptotics of the cross section when Q gets 
large but the dimensionless ratios xbj, z are held fixed, and q±/Q is either held fixed or 
integrated over. This has the unfortunate implication that the factorization theorem fails to 
give a useful result for the distribution in q± when q± <C Q. However, this is exactly where 
the cross section is biggest. Indeed, in lowest order, the hard scattering function has a factor 
S^(q±), while in higher order the hard scattering function has terms that diverge like l/q± 
times logarithms when q± — > 0. 

A convenient way of formulating these statements is to state that the factorization the- 
orem eq. (12) must be interpreted in the sense of distribution theory [26]. The theorem 
applies when the cross section is integrated with a smooth test function <f>(xBj, z, q±/Q). 
In the derivation [15,16,22], one applies to the hard scattering the approximation that the 
transverse momentum of the incoming parton can be neglected with respect to the transverse 
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momentum generated in the scattering, and one also neglects the transverse momentum gen- 
erated in the fragmentation. (Note that to the extent that these transverse momenta are 
not negligible, but are of order Q, the errors in the approximations are compensated by a 
correct treatment of higher order corrections to the hard scattering.) 

3.3 Factorization with intrinsic transverse momentum but no polarization 

To gain information on the q± dependence at small q±, we must derive a more powerful 
theorem that involves 'intrinsic transverse momentum' in both the distribution and the 
fragmentation functions. Such a theorem was derived for the Drell-Yan process and for the 
two-particle-inclusive cross section in e + e~ annihilation [27,17]. A similar theorem should 
apply here. An obvious ansatz is 

E ' E °wki = ? / rff / f / A - / ^ E ' E *»Jk * J 

+ Y(x Bi ,Q,z,q±/Q). 

(13) 

In this formula a represents the short distance part of elastic lepton-quark scattering. It 
contains a delta-function for momentum conservation. The sum over a is over all flavors of 
quark and antiquark. 

The first term on the right of eq. (13) dominates when q± <C Q. The second term, Y, 
is a correction term that enables eq. (13) to reproduce the ordinary factorization theorem 
eq. (12) at large transverse momentum, just as in the Drell-Yan case [17]. The Y term has the 
general form of the basic factorization theorem eq. (12), except that the low qj_ asymptote 
is subtracted from the hard scattering function. 

The function f a /A defined earlier gives the intrinsic transverse momentum dependence of 
partons in the initial state hadron. Similarly, Ds/ a gives the distribution of hadrons in a 
parton, with kb± being the transverse momentum of the parton relative to the hadron. 

Just as in [27,17], the hard scattering factor in the first term in eq. (13) can only be a 
2^2 process. Hence the fractional momenta of the incoming quark from hadron A and of 
hadron B in the outgoing quark are forced to be XBj and z. 

After integrating out the delta-function in a we obtain 

E ' EB wk~ B =^?E/^AmK,^) ^D B/a (z,k a± + q ± ) 
+ Y(x B j,Q,z,q±/Q). 

The picture that goes with these results is Fig. 2. All that we have done is to take 
account of the transverse momentum of the quarks relative to the measured initial- and 
final-state hadrons. This intrinsic transverse momentum has the effect of smearing out the 
delta function of q± that we remarked on earlier. The only generalization needed compared 
with the parton model is that the hard scattering can contain higher order virtual corrections. 

10 



(14) 



In the absence of gauge bosons in the strong interactions, this formula in the exact form given 
in eq. (13) is a theorem, that can be proved as in the Drell-Yan case [19]. 

The spin-1 gluons of QCD modify the theorem, by causing Sudakov form-factor effects. 
We expect that a proof can be given just as in the Drell-Yan case [17]. The effect is to 
broaden the transverse momentum distribution as Q increases, but in a spin-independent 
way: the broadening is due to recoil against the transverse momentum of soft gluon emission. 
This will have the effect of diluting the spin asymmetry we will discuss next. 



3.4 Factorization with intrinsic transverse momentum and polarization 



We now explain factorization for the semi-inclusive deep inelastic cross section when the 
incoming hadron A is transversely polarized but the lepton remains unpolarized. (It is left 
as an exercise to treat the most general case.) The factorization theorems, eq. (12) and 
eq. (14), continue to apply when we include polarization for the incoming hadron, but with 
the insertion of helicity density matrices for in and out quarks; this is a simple generalization 
of the results in [10,22]. 

The cross section will be linear in the transversity s± of the hadron (and also linear 
in its helicity A). Because transverse spin for a spin-! particle corresponds to off diagonal 
terms in the helicity density matrix, the other primary constraint comes from quark helicity 
conservation in the hard scattering, and this simplifies the factorization theorem. 

First, it is well known that at large transverse momentum, the transverse spin asymmetry 
is higher twist, as I now review. In that region, we use distribution and fragmentation func- 
tions integrated over intrinsic transverse momentum. Now, in the absence of a measurement 
of the polarization of the outgoing hadron, the single-particle fragmentation is spin inde- 
pendent. On the other hand, the transverse-spin dependence of the distribution functions 
is only in the off-diagonal elements of quark density matrices [8]. Therefore we need the 
part of the hard scattering that is off-diagonal in the helicity of the initial state quark but 
diagonal in the (summed) final-state helicities. Helicity conservation at the vertices for the 
gluon, photon and Z prohibits such a term, at leading twist. 

But, at low transverse momentum, the fragmentation function has dependence on trans- 
verse spin — see eq. (4). The corresponding hard scattering is just elastic electron-quark 
scattering, and we need terms that are off-diagonal in the final state quark helicity. Quark 
helicity conservation requires that these terms are also off-diagonal in initial-state quark 
helicity, and we are therefore discussing a spin transfer from the incoming to the outgoing 
quark. The off-diagonal terms in the density matrix for the incoming quark are given by the 
transverse spin distribution defined in eq. (2). 

Since there is a kinematic zero in the transverse spin dependence of the fragmentation at 
k± — 0, the spin dependence of the cross section must have transverse momentum dependence 
roughly proportional to the dimensionless coefficient 
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where M is a typical hadronic mass. This exhibits the kinematic zero when q± = 0, the 
leading twist asymmetry when q± = O(M), and the higher twist asymmetry when q± ^> M. 

We have now seen that the only transverse spin asymmetry at leading twist is in the low 
transverse momentum term in the generalized eq. (14). Moreover (5) implies that the trans- 
verse spin asymmetry has a sin dependence on the azimuth of the transverse momentum, 
so that any kind of uniform azimuthal averaging will remove the asymmetry. 

Following Meng, Olness and Soper [24] and earlier work, we could decompose the cross 
section in terms of scalar structure functions. In the case of fully inclusive unpolarized deeply 
inelastic scattering, with photon exchange, there are just two structure functions, the well- 
known F\ and F^. But when we measure one of the particles, B, in the final state, there is 
an extra vector in the problem, so there are more structure functions, enumerated in [24]. 
When in addition we allow the incoming hadron to be polarized, there are extra structure 
functions, just as for the kinematically isomorphic Drell-Yan cross section [10,28]. We will 
choose instead to work directly with the cross section and its angular dependence. It would 
of course be of interest to perform a structure function analysis. 

Since at leading twist, there is only transverse spin dependence in the low transverse 
momentum term in eq. (14), it is only this term that we need to change, with the result 

E Eb WW b = ^^J d ka± fa/A{xBi > ka±) paa ' ~d^~ D w( z > ^ + ^ (16) 

+ Y(x B j,Q,z,q±/Q). 

Temporarily we have changed notation: a, a', (3 and (3' represent indices for the initial 
and final-state densities matrices, p aa i is the helicity density matrix for the incoming quark, 
while the indices (3(3' on the fragmentation function D represent the dependence of the 
fragmentation on the outgoing quark's density matrix. 

We now recast this formula in terms of transverse spin vectors by using parity invariance 
as well as helicity conservation. First, we have already seen that helicity conservation tells 
us that the only transverse spin dependence is in the part of the hard scattering that is off- 
diagonal in helicity, i.e., <7+_ ;+ _ and its hermitian conjugate o"_ +; _ + . Thus the off-diagonal 
terms in the final and initial state density matrices are proportional: p°^_ = Cp+-, where C 
depends on the kinematic variables of the hard scattering. Furthermore, the expression for 
the off-diagonal term in a density matrix in terms of the spin vector is p + _ = 5|sj_|e* x , apart 
from a possible phase and a possible sign error in the exponent. Here x is the azimuthal 
angle of the transversity about the momentum of the quark. Furthermore parity invariance 
implies that if the spin of the initial quark is perpendicular to the plane of the scattering, 
then so is the spin of the final-state quark. Hence the coefficient C is real, if we choose our 
conventions such that a real value of p + _ corresponds to a spin vector perpendicular to the 
plane of scattering. 

It follows that the transversities of the initial and final quarks are proportional. The 
direction of the transversity of the final quark is obtained by rotating it about the normal 
to the plane of the hard scattering. (Fig. 3.) 

Hence the spin dependence of low transverse momentum term in the cross section is given 
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Fig. 3. Relation between spin vectors for initial and final state quarks. 
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by the following factor: 

l c E a e 2 a fdZk a± A T fAD(s>) 
E a e 2 Jd^k a jD 

Here we have used the fact that the spin transfer coefficient C, which we will calculate in the 
next section, depends only on the scattering angle of the hard scattering. It is independent 
of quark flavor and the other kinematic variables; thus we can factor it out of the integrals. 
The sums are over flavors of quark and antiquark. The transversity vector s' that is used 
in the spin-dependent fragmentation function AD(s') is the transverse spin vector of the 
initial hadron rotated in the center of mass of the hard scattering, according to the method 
of Fig. 3. 

The initial state distributions are functions of the length of k a ± only. So after integration 
over k a ±, the cross section acquires a spin dependence whose dependence on the azimuth 
of B is proportional to sin0, where is the azimuthal angle of ps, when measured in the 
lepton-parton center-of-mass about the outgoing jet axis (as defined by pj in eq. (10)), with 
= being the direction of the vector s' defined just above. 

3.5 Calculation of spin transfer 

We can easily compute the spin transfer part of the hard scattering cross section from the 
amplitudes given by Gastmans and Wu [29]. The cross section is for elastic electron-quark 
scattering, with all masses set to zero. Because of helicity conservation and parity invariance, 
there are two independent helicity amplitudes: M ++ and M + _, where the first and second 

indices label the helicities of the electron and quark. The amplitudes M and M_ + are 

related to these by parity invariance. 

The unpolarized cross section is a standard kinematic factor times | [|M ++ | 2 + |M + _| 2 ]. 
When the electron is unpolarized, the density matrix pf of the outgoing quark is expressed 
in terms of the density matrix p % of the incoming quark as follows: 

<>' - ( "r "r ) = (& 7 + ~ ) • 



where 



P-+ P- 



_ M ++ M* + _ + M_ + M*__ 

~ 12 i I M. 12 



M ++ | 2 + |M + _| 2 ' ^ 

Up to a phase convention, the explicit values for the matrix elements are [29]: 

M ++ = M__ = Aie 2 q a - 

1 — cos 9 . . 

(20) 

„ 1 -I- ms 

M_ + = M+_ = 2ie 2 q a - 



1 — cos 9 



Here q a is the electric charge of quark a in units of the basic charge e, while 9 is the scattering 
angle in the lepton-quark center-of-mass. The definitions of the helicity states are such that 
the scattering amplitude is purely imaginary; this convention differs from that of [29]. 
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It follows that the spin transfer coefficient is 



4(1 + cos 6) 

O — 



4+ (1 + cos#) 2 

(21) 

_ 1-2/ 



2 ' 



where y is the usual variable p A • (I — I') /p A ■ /, which is (E — E') /Em the target's rest frame 
and (1 — cos6>)/2 in the lepton-quark center-of-mass. The spin transfer coefficient is 100% 
at small scattering angles, and decreases to zero for exactly backward scattering. But even 
at 90° (y = 1/2), it is 80%. 

Note that the spin averaged square of the matrix element is 



4eX (, , 1 



M\* = —^(\-y + -y*\, (22) 



while the spin-dependent part is 



4eV 

C|M|2 = — ik(l-y). (23) 

y 



In the structure functions, the coefficient of 1 — y feeds into F 2) while the coefficient of y 2 /2 
feeds into F ± . Thus the transverse spin asymmetry we are discussing can be regarded as a 
100% asymmetry in F 2 and a zero asymmetry in F\. 



4 Predictions for DIS 

The measurements are of the single hadron distribution in the 'current quark fragmentation 
region' of collisions of wnpolarized leptons on transversely polarized nucleons. An asymmetry 
of the cross section is expected under reversal of the nucleon polarization as a function of 
the azimuthal angle of the outgoing hadron around the jet axis, which is defined by parton 
model kinematics. The asymmetry will have a sin <fi dependence. 

To a first approximation, the asymmetry is predicted to be the product of the polarization 
of the quarks in the initial nucleon, the spin transfer coefficient of the hard scattering, and 
the analyzing power of the fragmentation. In principle the asymmetry is a function of all 
the kinematic variables in the problem. The effect should be largest in the valence region 
for both the parton density and the fragmentation, say for x-q] ^ 0.3 and z > 0.3. One 
could even choose the hadron to be the leading hadron in the jet. 

The data should be analyzed separately for different charges, and ideally for different 
flavors of hadron (e.g., kaon versus pion versus proton). For scattering on a proton target, 
the biggest asymmetry will probably be for 7r + , since the 7r + has a valence u quark in common 
with the proton and deep inelastic scattering dominantly goes with u quarks in the valence 
region. It is possible that the asymmetry will reverse sign for tt~, and the value for n° is not 
obvious. It would also be interesting to see how big the asymmetry is for K + , which has a 
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valence u quark; this would probe the difference in the creation of strange quark- ant iquark 
pairs in fragmentation as compared to creation of u and d pairs. 

Furthermore there will be a strong dependence of the asymmetry on the transverse mo- 
mentum k± of the hadron relative to the jet axis. This dependence will be roughly of the 
form of (15), up to an overall coefficient, and it will be important to verify this dependence. 
At small k±, there is a kinematic zero, while, at large k±, QCD predicts the asymmetry to 
be of higher twist. It is at k±s of a few hundred MeV to about a GeV that one can expect 
the biggest asymmetry. 

Note that the unpolarized cross section itself has a significant azimuthal dependence [30]. 
This is predicted [31] to be leading twist at large k±, but higher twist at low k±, unlike the 
behavior of the transverse spin asymmetry. But clearly this would be a confounding effect 
to be careful about in the experiments. One must measure the actual spin asymmetry in the 
cross section and not just the azimuthal dependence of the cross section. 

Precise quantitative predictions cannot be made at present since the asymmetry depends 
on some not-yet-measured nonperturbative functions. Interesting comparisons can be made 
with data from proton-proton scattering and from final state correlations in e + e~ annihila- 
tion. Ultimately we must wait for experiments to give us the values of the nonperturbative 
functions. They will probe the spin and chiral structure of both the initial hadron and of 
the fragmentation. The spin and chiral structure of fragmentation is at present a very little 
explored subject. 

But we can suggest the general size, based on experience in other situations. The quarks 
in the hadron might be 50% polarized, as for the helicity asymmetry. The spin transfer in 
the hard scattering is mostly 80% to 100%, from calculation. Finally, the fragmentation 
could have an analyzing power of tens of percent. (The last figure is the most uncertain!). 
Overall we can then reasonably expect about a 10% asymmetry in the cross section. 

5 Hadron-Hadron Scattering 

Exactly analogous measurements can be made in jet production in hadron-hadron scattering. 
To define the jet axis, one could measure the jet axis in a conventional manner; this would 
probably have too large an uncertainty. An example of a better way would be to measure 
the correlation between the leading particles in opposite jets. Intrinsic transverse momen- 
tum effects generate out-of-plane transverse momentum for these particles. The maximum 
asymmetry will be when the spin vector of the incoming polarized proton is in the plane of 
the hard scattering. 

More theoretical analysis is needed here. 

6 Model Calculations of Fragmentation 

In this section I will present a very simple model calculation of the transversity dependent 
fragmentation function. It will show that the spin-dependence is permitted by the symme- 
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Fig. 4. Lowest order graph for fragmentation of quark to pion. The solid line is a quark and 
the dashed line is a pion. 

tries of QCD, but that it requires non-trivial phases in the interaction. Furthermore, it will 
be clear that the transversity dependence is caused by the spontaneous breaking of chiral 
symmetry. 

The model is a sigma model for pions coupled to quarks in the manner suggested by 
Georgi and Manohar [32]. The Lagrangian is to be considered an effective Lagrangian for 
suitable areas of nonperturbative QCD. While the model is not perfect, it does contain two 
important features of QCD: the quark degrees of freedom and the chiral symmetry breaking. 
We will calculate the lowest order graph for the fragmentation of a quark to a pion. To get a 
transversity dependence we will see that we must dress the propagators so that they acquire 
imaginary parts. 

The graph is shown is Fig. 4, and its contribution to the unpolarized fragmentation 
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follows from the definition eq. (3): 



u = ^ I dk ~j;dw) M+) ^ ~ pr ~ M2 ) (24) 

i tr 7 + (^ + M) 75 (^ -i> + M) 75 (^ + M). 

Here g is the pion quark coupling, presumably around a third of g^NN, for which 
9kNn/ '(16vt 2 ) ~ 1.1. The quark mass M is due to chiral symmetry breaking, so it should be 
around 300 MeV. For simplicity, we will neglect the pion mass. Then a simple calculation 
gives 

D = 16tt3 kl + M^ (25) 

and results in a transverse momentum distribution with the canonical 300 MeV width. 
The l/k\ tail at large kj_ is not to be trusted: it comes from a region where the quark 
is far off its mass shell, and there the sigma model is entirely inappropriate. Note the 
flavor structure: for different charge states in the fragmentation of a -u-quark, we have 
u — > 7r + : u — > 7T° : u — > 7T~ = 2 : 1 : 0. 



Similarly, the transversity dependent fragmentation function (eq. (4)) in our model is 
9 2 /■„- 1 

(26) 



i tr 7 + 75S± • 7± (A# + 5M) 75 (^ - + M) lb (A*$ + B*M). 

Here, we have anticipated the need for imaginary parts in the amplitudes, and have used 
a dressed quark propagator i(A]jt + BM)/(k 2 — M 2 ), where the scalar coefficients A and B 
depend on k 2 . On-shell renormalization gives them the value 1 when k 2 = M 2 . The result is 



g 2 2MlmA*B(l - z) 
16tt 3 {k\ + M 2 ) 
Hence the analyzing power is the ratio: 



^ = Ttz-* n„ 2 ^ M2 J v W - sW). (27) 



AD 2MlmA*B 1 - z ), „ IT T , lA 

= i ^( s » fc s - s x ± k y ). 28 

D /c 2 + M 2 z x 7 v 7 

This exhibits the properties claimed earlier: It gives a sin0 dependence on the azimuth 
of the transverse momentum, with a kinematic zero at k± = 0. A significant analyzing 
power evidently depends on the phases of A and B being substantially different, but there 
is presumably no reason why this should not be the case. 

There is an apparent blow up of (28) when z — > 0. This is not actually the case, since 
the virtual quark then approaches its mass shell, where A = B = 1. 
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7 Conclusions 

I have proposed that one can probe the transverse polarization of a quark that initiates a 
jet by measuring the azimuthal dependence of the hadrons in the jet, especially the leading 
hadrons. There are certainly other processes than the two explained in this paper that could 
be exploited. The typically large spin transfer in the hard scattering will aid in generating 
large asymmetries. 

Perhaps the most interesting aspect of the idea is that it gives one a direct handle on 
how the chiral properties of the strong interactions couple to fragmentation. In particular, it 
would appear from the model calculations in Sec. 6 that the fragmentation asymmetry will 
only appear if there is breaking of chiral symmetry. Thus, in the context of fragmentation, 
it probes the known spontaneous breaking of chiral symmetry. 

Conceivably, the measured asymmetry could be very small. This would require that the 
quarks in a transversely polarized proton have a small polarization or that the analyzing 
power of the fragmentation is small. The first possibility seems highly unlikely, given what 
we know of the other parton distributions in the valence region. The model calculation 
for the fragmentation indicates that there is no symmetry reason why the analyzing power 
should be small. Perhaps the best suggestion that a significant analyzing power is likely 
comes from the frequency of large spin asymmetries that are experimentally measured in 
hadron-induced processes. 

The differences between the longitudinal and the transverse spin asymmetries in the 
distribution function are sensitive to relativistic and spin-orbit effects in the proton wave 
function, and are thus of intrinsic interest. 
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Sivers [21] suggested that the k± distribution of a quark in a hadron could have an azimuthal 
asymmetry when the initial hadron has transverse polarization. However, as we will now 
show, such an asymmetry is prohibited because QCD is time-reversal invariant. 



Appendix 



Consider 



fa/A(x, k ± ;a,a r ) 




19 



which is a matrix in the helicity indices a and a'. This quantity is the same as eq. (1), 
except that we have replaced the hadron states by particular spin states in a helicity basis. 
(a and a' take the values L or R.) Contracting this matrix with the spin density matrix for 
the incoming hadron gives eq. (1). 

It can readily be checked that the matrix f a /A(x, k±; a, a') is hermitian, i.e., 

f a /A(x, k±; a, a') J = [f a /A(x, kj_; «)J = fa/A(x, k±; a, a'). 

Furthermore, let us choose coordinates in which the incoming hadron has zero transverse 
momentum, p± = 0. Then invariance under rotations about the z axis shows that the 
diagonal elements, f a /A{x, k±; RR) and f a /A(x,k±; LL) are functions of the length but not 
the angle of k±. Parity invariance implies that these two diagonal elements are equal: 

fa/A(x, k±; RR) = f a /A(x, k±; LL). 

We now apply a time-reversal transformation followed by a parity inversion. With the 
Bjorken and Drell [33] conventions for the Dirac matrices, we have 

{VT)^ a {x)VT = PTil>(-x), 

with PT = ?7°7 1 7 3 . Note that momenta are unchanged but helicity states get reversed under 
VT. Moreover there is a relative sign between the transformation of left and right-handed 
helicities: 

VT\p,L) = \p,R), 
VT\p,R) = -\p,L), 

where we have ignored an overall phase. The relative minus sign is essential, in order that 
T 2 be — 1 when acting on fermionic states. This sign can be verified from the transformation 
law of the Dirac field provided that one remembers that T is antilinear and that T* = — T. 

It follows that 

}2„ 



f a/A (x,k ± ;LR) = -J d y~^ e -^P + y-^ ± . y± 



(p, R\ T^ a (0, y~, y ± )T ^- T^ a (0)T \p, L)* 

djT dV ixp+y - +ik± . yj 
(2tt) 3 

(p, R\ MO, -y~, -y ± ) PT^-PTMO) \p, L)* 



~ -I 



= -fa/A(x,k±;LR), 

where the complex conjugation arises because the time reversal operator T is antilinear, and 
the last line follows from the hermiticity of /(««'). It is now immediate that the off-diagonal 
elements are zero: 

f(LR) = f(RL) = 0. 
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The matrix /(««') is therefore proportional to the unit matrix, so that there is no depen- 
dence of the transverse momentum distribution on the spin of the incoming hadron. This 
contradicts Sivers' [Sivers] suggestion. 
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